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Abstract 

An iterative scheme is set up for solving the loop equation of the her- 
mitian one-matrix model with a multi-cut structure. Explicit results are 
presented for genus one for an arbitrary but finite number of cuts. Due 
to the complicated form of the boundary conditions, the loop correlators 
now contain elliptic integrals. This demonstrates the existence of new uni- 
versality classes for the hermitian matrix model. The two-cut solution is 
investigated in more detail, including the double-scaling limit. It is shown, 
that in special cases it differs from the known continuum solution with one 
cut. 



1 Introduction 



The study of multi-cut solutions of matrix models is interesting for sev- 
eral reasons. First of all it completes the perturbative analysis around 
N = oo, exhausting the full parameter space of the coupling constants for all 
possible kinds of solutions. This is important in order to tell, whether there 
exist different types of phase transitions, universality classes or continuum 
limits than the ones, which are already known for the one-cut solution of 
the model. 

Secondly there may be applications apart from quantum gravity to two 
dimensional QCD or statistical physics. Results of random matrix theory 
are used for example in the study of the spectrum of the QCD inspired dis- 
cretised Dirac operator. The merging of several cuts into a one-cut solution 
of a complex matrix model with "temperature" would be of interest there 
in various limits (|l| and references therein). 

In this paper the approach of iteratively solving the loop equation for 
the hermitian one-matrix model [Q] is chosen and generalised to an arbi- 
trary but finite number of cuts. The success of this method is based on 
a nonperturbative treatment in the coupling constants, which is valid for 
an arbitrary polynomial potential. Even away from the double scaling limit 
higher orders in the genus expansion can be obtained explicitly. The validity 
for an arbitrary potential makes it possible to read off immediately, whether 
the results are universal or not. Here this will lead to a whole set of new 
universality classes for the planar two-point function or two-loop correlator. 

Another reason for choosing the loop equation techniques is that it seems, 
that the method of orthogonal polynomials is not applicable to multi-cut 
solutions in general. Solving the string equation for the norms of the poly- 
nomials numerically, instabilities or a so-called chaotic behaviour has been 
found by a variety of authors []|, ||], ||, ||, [j]]. The region, where these instabil- 
ities occur, seems to match precisely with the part of the phase space, where 
the multi-cut solutions dominate. This has been shown for the hermitian 
model with an even sextic potential as an example @, [||]. So apart from 
the symmetric two-cut solution and some cases of degenerate minima of the 
potential, where smooth numerical solutions have been found |3], ||, ||], an 
analytic expression for the orthogonal polynomials even for the planar limit 
is still missing in the case of a generic multi-cut situation^ As an example 

1 The two- component ansatz for the recursion coefficients made in 0], |n| and is 
restricted to the purely symmetric case of two cuts. The failure of a more general ansatz 
for more cuts is studied in lid]. 
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the nonsymmetric two-cut solution, which is presented here most explicitly, 
explores such a region of chaotic numerical solutions 

The paper is divided up into two parts. In the first part the iterative 
scheme for solving the model with s cuts is presented, including results for 
genus one. The second part is then devoted to a more detailed study of the 
two-cut solution, which also contains an analysis of the double-scaling limit. 

In section 2 and 3 the main definitions and the loop equations for the 
hermitian model with multiple cuts are displayed, following closely the nota- 
tion of Q. The planar solution is derived in section 4, including the required 
boundary conditions. Then, in section 5 the iterative scheme for higher gen- 
era is presented. Special care has to be taken to the inversion of the loop 
equation due to the zero modes, which occur for more than one cut. The 
one-loop correlator of genus one is given for an arbitrary number of cuts. 

Turning to the two-cut solution, section 6 contains the full two-point 
function as well as the free energy of genus one. After a short discussion 
of general properties in the symmetric limit, section 7 deals with the differ- 
ent possibilities of taking a continuum limit, displaying some new results. 
Section 8 closes with a final discussion and future prospects. 



2 Basic definitions 

Throughout the paper the same notation as in ^ is used, which is re- 
displayed here for completeness. The partition function of the hermitian 
one-matrix model is defined by 

Z[N,{ 9i }] = e^W*}] = J d<\> exp(-iV TrV^)) , (2.1) 

where the integration is over hermitian N x N matrices <j). The matrix 
potential is given by the following power series, 

oo 

V{<t>) = E-^' • (2-2) 

A specific potential of finite order can be inspected by setting the extra 
couplings to zero in the final result. In this way of keeping all the coupling 
constants gi throughout the calculation, they can be used as sources for 
expectation values^] of the following type, 

- m J-F [N,{ 9i }} = -J-<Tr<F> , mGN + . (2.3) 
dgm N 

2 Averages are defined as usual by {Q((f>)) = ^ J d<f> Q(4>) exp(— JV TrV(^)) . 
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Introducing the loop insertion operator 

the generating functional or one- loop average W(p) can thus be obtained 
from the free energy F 

{ ' ~ N^ p k + l N\ p-<f> 

= ^(p)F+-. (2.5) 
av p 



In the same way all the multi-loop correlators can be derived by applying 
$(p) to F (or to W(p)) 

1 „ 1 



W( Pu ...,p n ) = N n ~ 2 ( Tr — —-j • • • Tr 



Pl~<t> Pn-<t> 

= ^(Pn)^(Pn-l)---^(Pl)F, n>2. (2.6) 

Here conn refers to the connected part. As the loop correlators and the free 
energy have the same genus expansion, 

F = Eiv^' (2 - 7) 

9=0 

oo 

W{ Pl ,..., Pn ) = ]r_W 9 ( Pl ,..., Pn ) , (2.8) 

3=0 

eq. (^) is valid for each genus g > separately. Finally the asymptotic 
behaviour of W(p) for large p can be read off from the definition ( [2.5| ) 

lim W(p) ~ - . (2.9) 

p— >±oo p 

The r.h.s. does not depend on N, so the leading contribution clearly comes 
from the planar part 

lim W (p) ~ - , 

p^±oo p 

lim W ff (p) ~ O(l) , $ > 1 . (2.10) 
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Using this fact the genus expansion of eq. ( |2.5[ ) reads 
W g (p) = 4?(P) F 9 > 9 > 1 , 



(2.11) 



so the p-dependence of W g (p) is completely absorbed in the total derivative. 
The last two equations will become important, when the genus expanded 
loop equation is inverted, which determines W(p) iteratively in genus. 

3 The loop equation 

The derivation of the loop equation for multiple cuts can be performed 
exactly along the same lines like for the one-cut case, exploiting the invari- 
ance of the partition function under a field redefinition <p — > <f> -\- e/(p — (j)). 
The result differs only by the contour C of the complex integral (s. fig. [l]) 

£ &h w ^ = (M '» )2 + ww^ w <>> <**'< (3 - 1) 

where V'(uj) = J2j ■ I n the derivation it has been assumed like in the 

one-cut case (e.g. in [jl3|), that the density pjv(A) = <5(A — At)) of 

the eigenvalues of the matrices (j) has a compact support a in the vicinity 
of N = oo. Here a consists of an arbitrary but fixed number s of distinct 
intervals 

s 

a = {J[%2i,X2i-l] , x 1 > x 2 > ■ ■ ■ > x 2s . (3.2) 
1=1 

From rewriting W(p) in terms of PnW 

W(p) = I dX^l , (3.3) 
J p - A 

it is clear, that when expanded at N = oo W{p) has s cuts along the real axis 
on a and is analytic elsewhere. Hence the contour C in the loop equation 
( |3.1| ) has to enclose all singularities of W(co), but not the point u> = p. It 
may now well be situated between two cuts, in particular a double-scaling 
limit can be performed at an internal edge of a cut (s. chap. |7|). 

Inserting the genus expansion eq. ( [2.8D into the loop equation (|3.l|) one 
gets 

r ^v^) = 2 

Jc 2ni p — lo 
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Figure 1: The contour of integration C = U| =1 Cj 
for genus zero and 

g-1 , 

(IC-2Wo(p))W g (p) = J2W gl (p)W g - gl (p) + —(p)W g ^( P ), g>l, 

g/=l " 

(3-5) 

for higher genera, where AC is a linear integral operator defined by 

tm s I ■ (3.6) 

Jc 2iri p — ijj 

Knowing the result for Wq{ P ) one can determine W g (p) for g > 1 iteratively 
from contributions of lower genus on the r.h.s. of eq. ([3.5]), provided that 
the operator (K — 2Wq(p)) can be inverted uniquely. 

4 The planar solution 

The quadratic equation ( [Q| ) for the planar solution Wq( P ) can be solved 
in the following way. Deforming the contour in eq. ( |3.4| ) to infinity one gets 

(W (p)) 2 = V'(p)W (p) + i P:^-W (u). (4.1) 

J Cao 2mp-uj 

The solution is then formally given by 

w (p) = \v'{p)-\^{v'( P )y + ±Q{p) , 
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Q(P) = I ^^oM, (4.2) 
J Coo 2-Kip-uj 

where the minus sign is chosen in order to fulfill the asymptotic eq. fl2.10|) . 
Now if it is assumed that Wo(p) has s cuts in the complex plane, one can 
make the following ansatz for the square root, 



W (p) = \lv'(p)-M(p) ] IY[ 2 " =i (p-x l )J . (4.3) 

M{p) is assumed to be an analytic function, which still has to be determined. 
The signs of the complex square root are defined by the following choice, 



r 2s 

\i I 

p— >±oo 



inn , [J. ^(p - Xi ) ~ p s . (4.4) 



The analyticity of M(p) allows one to write it as 



J Coo 2muj-p 



Now solving eq. ( [4.3| ) for M(u) and plugging it into eq. ( f4.5| ) yields 

M(p) = f — , \ ! , (4.6) 

because the Wo(u))-term vanishes due to its asymptotic behaviour at infinity 
~ ~. When evaluated for a given potential of finite degree d, M(p) will be 
a polynomial of degree d — s — 1. Therefore the maximal possible number 
of cuts has to be s max < d — 1 0. Reinserting the solution for M(p) into eq. 
( [4.3D and deforming back the contour leads to the final result, the planar 
one-loop correlator with an s-cut structure, 



1 l dw V'{uj) 

W (p) 



2 Jc 2m p — lo\^ 



n(^)- («> 



From this result the planar eigenvalue density p(A) = liniAr^oo pat (A) can be 



recovered, taking the limit N — > oo of eq. (3.3) 



p(A) = -^Hm(Wo(A-^)-W (A + ^) 

2m e^O V 



i-|M(A)|^n^(A^) , AG a. (4.8) 



3 It can be shown, that even s m ax < i holds. 
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Equivalently p(X) could have been obtained from the saddle point equation 
for the partition function in terms of the eigenvalues (see for example in 
|p.2|). Here the solution takes a closed form for an arbitrary number of cuts 
s. 

In the planar solution for Wq(p) eq. ( |OD the edges of the cuts 
i = l,...,2s , still have to be determined as functions of the coupling 
constants gt, i G N + . These boundary conditions are derived by exploiting 
again the asymptotic behaviour of Wo(p) ~ h- All the coefficients in eq. 
(|4.7|) of the order p s_1 down to p° must vanish, whereas the ^-coefficient 
must equal unity, 

. 1 f du uj k V'{uj) . n n . 



2 Jc 2m 



These conditions provide only s + 1 equations for the 2s constants Xi . This 
is sufficient if and only if s = 1. For s > 2 the remaining equations are given 
by a criterion of stability against tunnelling of eigenvalues between different 
cuts in the limit iV — > oo [p^| . Stability is achieved by imposing equality 
of the chemical potentials of all cuts in the saddle point equation for the 
density p(\). These conditions can be most suitably written in the form [|l5|] 

%2k 



0= j d\ M{X)^l[^ =i (X-x t ) , k = l,...,s-l, (4.10) 



x 2k + l 



where the integration is performed between each two neighbouring cuts. 
This second type of boundary condition for s > 2 leads to the appearance 
of elliptic integrals in the higher genus correlators as well as in the planar 
multi-loop correlators. 



5 The iterative procedure 

After having solved the planar part of W(p) as a starting point, higher 
genus contributions can now be iteratively determined by inverting the genus 
expanded loop equation ( |3.5| ). From W g (p) all multi-loop correlators of the 
same genus can be obtained then by simply applying the loop insertion 
operator -jw(p) to it (see eq. ( ^lf ) and (2.5)). Along the same lines like in the 
one-cut solution B a change of variables from coupling constants to moments 

(k) 

Mf' will allow to express the higher genus correlators nonperturbatively in 
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the coupling constants gi, depending then only on a finite number of these 
moments, 



Afi fc) 



f^V'i^icj), kN + , i = l,...,2s, (5.1) 



with 



f ( W ) ^ r -L 13f ^(c) , 0W ( w ) = . (5.2) 



Deforming the contour to infinity and expanding $\ k \uj) the moments can 
be seen to depend in the following way on the couplings, 



2s , 

M i k) = gk+s + gk+ s +iJ2xj{- + kdij) + ... . (5.3) 

3=1 

Here an explicit dependence on the number of cuts s enters. The moments 
may also be used to characterise multi-critical points. Namely because of 

the m-th multi-critical point at x% is given by 

M i (fe) = 0, k = 1, ...,m-l , 

M/ m) ^ 0, (5.5) 



which indicates that the eigenvalue density p(\) in eq. ( |4.8D accumulates 
m — 1 extra zeros at Xj. 

Secondly it can be seen from eqs. ( |5.4[) and ([4.3]), that the planar one- 
loop correlator depends on all moments, when M(A) is expanded in a Taylor 
series. This is a special feature of genus g = only. 

5.1 Determination of the basis 

In order to proceed it is necessary to introduce a basis for the operator 



acting on W g {p) in eq. (3.5) 



{t - 2W (p)) xS n) (p) = - , n € N+, i = 1,... ,2s . (5.6) 

{P Xi) 
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Given that the r.h.s. of eq. ( |3.5| ) is a fractional rational function of p having 
poles at the Xi only, W g (p) will then have the following structure 

3g-l 2s 

W 9 (p) = EE^Sx^b), g>l. (5.7) 

n=l i=l 

(n) (k) 

The A\ J are complicated functions of the Xi and the moments M- . As the 
order of the highest pole in W g (p) is not changed by assuming a multi-cut 
structure, W g (p) will depend on at most 2s(3g — 1) moments, following the 
same arguments like in the one-cut solution [g]. 

A set of basis functions fulfilling eq. ( |5.6| ) is defined by 

xt\p) E 

Mf \ k=l / 

(5- 

which can be proven by induction. However, this definition is not unique, 
as the kernel of (/C — 2Wo (p) ) is not empty, 

Ker(/C - 2W (p)) = Span{j?V (0) (p); I = 0, . . . ,s} . (5.9) 

This can be shown by using the boundary conditions of the first type eq. 
fl4.9| ). Now because of the asymptotic behaviour eq. ( |2.10| ) only terms 
~ O(^z) at large p may be added to W g (p) for g > 1. This requirement 
reduces the set of zero modes to be 

allowed zero modes : {p l 4>^°\p)] I = 0, . . . , s — 2} . (5.10) 

In particular for the one-cut solution this argument excluded any zero mode 
to be added, making the definition of the basis in eq. ( |5.8| ) unique. In the 
general case any linear combination of the s — 1 functions can in principle 
be added to W g (p) for a solution with s cuts. 

In addition to its asymptotic W g {p) has to fulfill eq. ( [2.11 ), which can be 



used to fix the basis uniquely. It says that the p-dependence of W g {p) must 
be completely absorbable into derivatives with respect to -^{p)- This will 
only be possible, if the basis functions xt^ (?) ma y be expressed completely 

in terms of ^t-(p) and d y (p) as functions of p. So the x\ (?) m e Q- ( |5.8| ) 
must be redefined in order to achieve this, which will precisely fix its zero 

(k) 

mode content. The derivatives of the Xi and the moments M- can be 
obtained by applying the loop insertion operator -^(p) to the definitions 
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( |5.1| ) and to the boundary conditions, when rewriting it in the following 

way: 

d d 2s dx- d 2s 00 dM (fc) d 

i=l 1 i=l k=l om i 

f) CO ■ o 

Using the identity 

the result for the moments reads 



(P) = (^ + ^)(Mf +1) ^(p)-# +1) (p)) 

+ 2 ^ f- fx, - a^fc-i+i w {P, J 



dV 

2s k 



*U&w- x ir- t+1K 1 dv 

* = !,..., 2s , G N + . (5.13) 



The quantities ^ (p) are given by the solution of the following set of linear 
equations, 

= jr^M^^ip) - +2*p fc -V ) (p), fe = o,..., a> 

= E (V^) - ^(p)) *u > j = 1, . . . , s — 1 , 

^, = y dA (A _ x0 , (5.14) 

where -^ip) has been applied to eqs. ( fO| ) and ( [4 . 1 0| ) . The result of the 
latter is derived in appendix [A], where some care has to be taken to the 
interchanging of -§j(p) and deforming the contour in eq. ( [4.6[ ) for M(X) to 
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infinity. It can easily be seen from the linear system of equations ( |5.14 ), 
that the solution will always take the form 

mP^(p) = <t> ( -\ P ) + E<W0 (O) (P) . i = l,...,2a . (5.15) 
av 1=0 

The aij, i = 1, . . . , 2s, I = 0, . . . , s — 2, will then be determined by the 
following set of equations: 

o = EXXizfp' - EeV 1_ V + 2k P k -\ k = o,..., s , 

4=1 (=0 i=l «=0 

2s s-2 

= EE^^iP 1 ' j = l,...,s-l. (5.16) 

4=1 1=0 

Comparing the coefficients in powersQof p , 1 = 0,...,, s—2, leads to 2s(s— 1) 
equations for the j. They will only depend on the X{ and K^j. Hence the 
only difference between the ^jy(p) in eq. (|5.15 ) and the respective quantities 
in the one-cut solution Q is the occurrence of a linear combination of the 
allowed zero modes from eq. ( |5.1C ) . More explicit results for s = 2 are given 
in the next chapter. 

Now all necessary ingredients for redefining the basis (|5.8|) in terms of 
total derivatives are collected. Solving eq. ( |5.13 ) for 0$ (p) and express- 
ing its p-dependence as ^(p)-terms plus zero modes, a unique basis can now 
be obtained inductively from eq. ( |5.8| ) by subtracting these corresponding 
zero modes: 



M> 



dv-P art 



k=l 



i = l,...,2s , n€N+ . (5.17) 
The first basis functions then read 

xf\p) = ^(p), i = l,...,2s, 

2 d 1 ^* d 

X ^ ip) = ~3dy (p)ln|Mi(1)| ~ 3^dy (i?)ln|xi " ■ (5 ' 18) 

3=1 



'The coefficient of p fc 1 is zero for k = s. 
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5.2 Calculation of genus 1 

Having determined the basis, the loop equation can now be inverted step 
by step in genus. For genus g = 1 eq. (^) reads 

(£-2Wb(p))Wi(p) = jL(p)W (p). (5.19) 

Using the result for Wo(p) eq. ( |4.7|) and the loop insertion operator from 
eq. ( 5.11 ) the r.h.s. is given by 

d 3 if, 1 1 2 " 



dy (p) Wb(p) 16 § (p - x,) 2 8 £ (P " *i)(P " *i) 



1 1 ^i, 1 (i) dxi 

J_ y, 1 1 y, 1 / _L 1 \ 

16 ^ (p - x;) 2 8 xj - y (p - Xi) (p - Xj) J 

-. 2s s-2 l 
1 «^ 

4^^ P-x, 

Here the fact has been used, that all the regular parts coming from p p _ x , 2 = 

1, . . . , s — 2, will vanish due to eq. ( |5.16| ), such that Wo(p,p) = ^(p)Wo(p) 
fulfils its correct asymptotic. The two-loop correlator at different arguments 
Wo (p, g) can be obtained in the same way, which is derived in appendix [l|. 
The result for the one-loop correlator of genus one with s cuts can now easily 
be obtained by using the basis eq. ( 5.18 ), 

i=l i,j=l 1 i 

+lEE«M^xS 1) (p) • (5.21) 

i=l 1=0 

The integration of W\(p) in order to get F± as well as the calculation of 
higher genera gets technically very much involved, as the integrals Ki j then 
have to be integrated or differentiated with respect to the Xj, i = 1, . . . , 2s. 

For the case of two cuts, the integrals K^i, i = 1, . . . , 4 , can be expressed 
by the well known complete elliptic integrals of the first, second and third 
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kind, making a more detailed analysis possible. This will be the subject of 
the next chapter. 



6 The two-cut solution 

The case of two cuts may already appear for the symmetric quartic po- 
tential, when the coupling constants are suitably chosen |l6| . The approach 
presented here allows for a closed treatment of an arbitrary polynomial po- 
tential including higher genus contributions. It is beyond the scope of the 
method of orthogonal polynomials, as it has been mentioned already in the 
introduction. 

The explicit solution for the zero mode coefficients ohq i n ec i- 

(|5l5l) will 



allow to inspect the planar two-loop correlator in more detail. Knowing the 
complete Xj-dependence of these coefficients makes it also possible to study 
the double-scaling limit. The set of equations ( |5.16 ) reads for s = 2 



i 



= 

i=l 

4 

= ^ctiXi - 2 
i=l 

4 

= J2 ( a i X i 
8=1 

4 

= Y. a i R ii 



8=1 



where on = a^o and 



h '> = I !lX " (\ - Xl ) ' i = 1 .---. 4 - ( 6 - 2 ) 



The solution for the integrals K{ in terms of complete elliptic integrals may 
be taken from |L7]], the result being displayed in appendix |C|. The solution 
of the eqs. ( |6.1D reads 

M^^ip) = 4\p) + a^ (0) (p), i = l,...,4, (6.3) 
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where the cti can be found after a tedious calculation to be 
ai = ( 1 + 

X\ — Xi V 

a 2 = f 1 + 

a 3 = ( 1 + 

X3-X2 \ 

4 X4 — X\ ( ^ (;(•;] - /V(/rj 

fc 2 = Ql -X A ){X 2 -g 3 ) /g^ 

(xi - x 3 )(x 2 -x 4 ) 

Having the explicit result for the ^p-(p) at hand, the planar two- loop cor- 
relator can be evaluated also for different arguments most explicitly. It is 
derived in appendix [B], which reads 



{X2 


- x 4 ) E(k) 


(xi 


- x 2 ) K{k) 


(xi 


- x 3 ) E(k) 


{X2 


- Xl )K{k) 


{X2 


- x A ) E{k) 


(x 4 


- x 3 ) K(k) 


(xi 


-x 3 )E(k) 



W (p,q) 



1 



4(p - qf 



+ 



l(p- 


xi)(p - 


x 4 )(q - 


X2){q - 


X3) 


(p - 


x 2 )(p - 


X3)(q - 


x\){q ~ 


Xi) 


Kp- 


- x 2 )(p- 


- x 3 )(q - 


- xi)(q - 


- X4) 



(p - xi)(p - Xi)(q - x 2 ){q - x 3 ) 



+ I / 4 1 " 9T~~^2 (6-5) 

4 Vnli(p-^)(9-^)^) 2 (^-^) 2 



The function /({xj}) of the X{ is given in eq. ( |B.6| ). The result fulfils its 
required analyticity properties (see appendix |B|) . As it must be regular at 
coinciding arguments, it has to be compared with eq. (|5.20 ) for the two-cut 
solution, which can also be found in appendix M, 



1 A 1 1 / 1 



m "' p) = u £ + 4\(p-zi)(p- ») + (p - i55 - ») 

lA 1 11 

(p-^b-x,) + 4n 4 =1 (p-x,)^) /(|Xj}) • (6 - 6) 

These results for the planar two-loop correlator are clearly universal, as they 
depend on the coupling constants only implicitly via the endpoints of the 
cuts. The part proportional to the elliptic integrals simplifies considerably 
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in the limit of a symmetric potential, where X4 = —x\ and X3 = — X2- The 
modulus will then be 

2 4x l3 ; 2 , . 

and the function /({xj}) becomes 

/ S j/m({^i}) = (xi + x 2 ) 2 . (6.8) 

6.1 Results for genus g = 1 

With the more explicit results for the two-cut solution at hand the free 
energy F\ can now be integrated from W\{p). As it has been mentioned 
already the knowledge of F\ in not necessary to get the higher orders of 
Wg{p) though. Inserting the basis functions eq. ( |5.18|) into eq. Q5.21 ) W\{p) 
reads 

1 4 d Id 

w i(p) = -^E^W^^VgE^WM^-^l 

i=l i<j 
1=1 

Using the well known relation for the integral K(k), 

^K(k) = ^h^(m - (1 - , (6.10) 

the following helpful relation can be derived, 

d ^ ctx ' d d 

2— (p)]n\K(k)\ = -^a i ^(p) + —(p)ln\x 1 -x 3 \ + — (p)ln\x2-x 4 \. 

(6.11) 



The free energy F\ now can be read off easily from eq. (2.11), W\(p) = 

i=l j<jr 

+ j^ln |xi - x 3 | + In |x 2 - a: 4 |^ . (6-12) 
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Taking the symmetric limitf], the result looks very similar to the free energy 
of genus one of the 0(n)-model for n = +2 

The calculation of higher genera may now be performed by computer 



algebra, using eq. (3.5). Taking the derivative jy(p), all quantities can 



expressed algebraically by the known expressions for the ^(p), d y (p) and 



K(k) and E(k), using eq. (6.10) and a similar expression for E(k). 



The particular case of having two cuts allows for a certain check of 
the one-loop correlator, which is valid for all genera. If one is taking the 
symmetric limit at the end of the calculation, which is implemented by 
setting the odd coupling constants g2i+\, i € N, to zero, all the odd powers 
of expectation values of matrices will vanish, 

(Tr0 2fe+1 ) S2/m = keN . (6.13) 

Consequently W(p) will become an odd function in p, 

1 1 (Tr<p 2 ) svm 
= - + — V , ym + ... . 6.15 

p TV p A 



In particular the term ~ O(^z) will disappear for all genera. Now the only 
function of p contained in the basis with this asymptotic behaviour^] is the 
zero mode (p^°\p). Hence it is not contributing to the symmetric one- loop 
correlator Wg ym (p) for all genera. Using this special property the result for 
W\(jp) can be checked in this limit. Extracting the 4>^ (p)-dependence from 
eqs. ( gj) and ( PH) it can be shown, that its coefficient indeed vanishes in 
the symmetric limit. 



3 The moments will then behave like M[ n) = (-)"M 1 (n) , M< n) = (-)™M 2 Cn) . 
3 For s > 3 there will always be zero modes left with an asymptotic ~ 0(\). 
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7 The double-scaling limit of the two-cut solution 



It has been argued that in the double-scaling limit (d.s.l.) all the multi- 
cut solutions should be equivalent to the one-cut solution for the following 
reason. Adjusting the coupling constants to achieve a multi-critical be- 
haviour at one of the edges of the cuts, say at Xj, the d.s.l. magnifies the 
accumulation of zeros at Xj such that the structure elsewhere can be ne- 
glected. However, this argument is only true, if the cuts do not touch or 
vanish at the same time. The given solution with two cuts therefore supplies 
an example, where this can be checked explicitly. It will be verified, that the 
generic d.s.l. at any Xj is equivalent to the d.s.l. of the one-cut solution ||. 
But if the scaling limit is taken at a cut, which simultaneously shrinks to 
zero, or merges with the other cut, or both, a different continuum behaviour 
will turn out. 



7.1 Scaling limit at Xj 

In the generic case the d.s.l. is performed at one specific edge of the cuts 
Xj, where for the m-th multi-critical point m — 1 extra zeros accumulate in 
the eigenvalue density eq. Q4.8P . Fixing the coupling constants appropriately 
to reach this point, p and Xj will scale in the following way 

Xj = Xj — aA™ 

p = Xj + an , (7-1) 

whereas the x^j are kept fixed. For an m-th multi-critical point the mo- 
ments scale according to 

Mf ] ~ a m - k , fc = l,...,m-l, (7.2) 

!k) 

where again the M-A do not scale. Looking at the explicit solution for the 



first basis function, which is given by eqs. fl6,3| ) and (|6.4|), it is clear, that the 
zero mode contribution will be sub-dominant in the scaling limit eq. (|7.1|) 



%ip) = T^^iP) (d-Bl) • (7-3) 

3 



The derivatives of the moments take the form 



dV ^' v 2' V 3 dV TJ 



(P) = ik + =-)(Mi k+1 ^(p)-^ +1 \p) 
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4 +1) (p) = 7 1 ]k+l Tl x 1 = (d.s.1.) • (7.4) 

Looking at the inductive construction of the basis eq. (|5,17|) it is obvious 
then, that the zero modes will be also sub-dominant for the rest of the basis, 
which is hence given by 

xf ] (p) - f4 n) (p)-E^ +1) xf (P)) . (d.s.1.) 











CL~ m ~ n 


4 



k=l 



(7.5) 



Consequently in the scaling limit eq. ( |7. 1| ) the basis for the poles at Xj 
precisely looks like the one for the one-cut case at x = Xj ||. The other 
basis functions, which are sub-dominant ~ az, will not be needed, as only 
contributions from poles at dominant. The complete equivalence to 

the one-cut solution in the d.s.1. will now be shown by looking at the starting 
point of the iteration Wo(p,p) and at the scaled loop insertion operator. 
It is easy to see from eq. fl5.20 ), that the leading part in orders of a will 



be 

w "^ = h<^r (d • s ■ L, • (7 - 6) 

which is the same like in the one-cut solution. The scaled -wj(p) from eq. 



(5.11) will take the form 



d , , dxj . . d ^ dM^ d . . . . 



dx dM {k) 



where -^t(p) and (p) are given by eq. ( |7.3| ) and eq. (7.4) respectively. 
Hence the generic d.s.1. of the two-cut solution precisely maps to the d.s.1. of 
the one-cut solution in ^ , when the quantities d c and M/% there are replaced 
by d c = ]J i -^j(Xj — Xi) and M- respectively here. However, one has to bare 
in mind, that in this symbolical equivalence of the respective formulas the 
moments will explicitly depend on the number of cuts as functions of the 



couplings (see eq. (|5.3|) ). 

In reference Q explicit results up to and including genus g = 4 are 
presented in the d.s.1., which are valid here as well. The analysis concerning 



the possible kinds of combinations of moments, that can appear in eq. (5.7) 
also applies here. 



18 



7.2 Merging and shrinking cuts in the d.s.l. 

A different scaling behaviour may be expected, when the d.s.l. is per- 
formed at a point, where the two cuts merge, or one of them shrinks to 
zero, or both happens together. When the case of merging is considered, it 
has to be taken into account, that extra zeros of M{p) lying between the 
cuts will be picked up in this limit 7 . This will yet only affect the order of 
multi-criticality, but not the scaling behaviour of the zero modes. 

The outcome of the analysis will be, that in the d.s.l. at either merging or 
shrinking cuts the zero mode contributions will stay sub-dominant, although 
becoming singular in the latter case. Nevertheless the continuum results will 
change in all cases, as in eq. ( 5.2C| ) more terms will contribute. When the 



shrinking and merging of the cuts is considered simultaneously, say X2, #4 — > 
X3, even the zero modes will be enhanced such that they contribute in the 
d.s.l., leading to another different continuum theory. 

First the situation of the two cuts merging is considered in the limit 
%2 — * %3, which can be parametrised by 

X2 = £3 + av , v > , 
x 3 = xl — afi, n > , 

p = xz + air , (7.8) 

as at the critical point p — x% will also scale. Looking at the explicit results 
eq. ( |6.4| ) for the zero modes, the modulus k 2 will be of the order a, 

So the elliptic integrals can be expanded in k 2 , which reads 

§<| = l-^ + CV). (T.10) 

The scaling behaviour in 02 and coming from the factor X J_ X3 will thus 
be cancelled, leaving them sub-dominant as in the generic case, 

1 



Oil 



X\ — X3 



7 From analyticity properties and the positivity of the eigenvalue density it follows H], 
that M(p) will always have at least one real zero between two adjacent cuts for arbitrary 
s. 



19 



n ' 2 ! 2{xl-x A ) 

1 (xi - x 4 ) 



x% 




X\ 




1 




X 3 




X,\ 




1 




:r 4 




x 3 



" * * 1 2(xi -xjj 

«4 = ^ • (7.11) 

X4 — X\ 

The coefficients ct% and Q4 are regular in the limit eq. ( [7.8[ ) anyway. The 
result for the ^p-(p), i = 2,3 , will therefore be the same like in eq. (|7.3|) . 
Still, the continuum limit is altered, as in eq. ( |5.13| ) more terms will survive 
now in comparison to eq. ( [7.4[) in the generic case, 

and the same for the indices 3 and 2 interchanged. The d.s.l. performed 
at the interior of the two cuts merging therefore leads to a new continuum 
limit of the theory. The starting point for the iterative solution in the scaling 
limit reads 

w (p,p) = ± £ 7-^2 - l—^— (— — ) • ( 7 - 13 ) 

16 -r^ 3 [p -Xi) z 8 x 2 - x 3 \p - x 2 p-x 3 J 



Using eqs. ( |7.3| ) and ( 7.12j ) to obtain the double-scaled basis, the genus one 



result can easily be achieved. 

Second, the case of the scaling limit is considered, which is performed 
at a cut simultaneously shrinking to zero, say X4 — ► X3. It is conveniently 
parametrised by 

x 3 = ^3 + av 1 v > > 
X4 = X3 — a\i , fj, > , 

p = xg + avr, (7.14) 

where the other possible case x 2 —> x\ can be obtained by interchanging the 
indices 4 <-> 1 and 3 <-> 2 everywhere. The modulus k 2 will reach unity in the 
limit ( [7.14 ), where the integral of the first kind K(k) becomes singular. The 
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integrals in eq. ([O]) can now be expanded in terms of the complementary 
modulus A/ 2 , 

u /2 _ -, , 2 _ „/„ , „\ ( ; ' \ rnf„2\ 



k = 1 — k = a(u + M ) + Q(a 2 ) . (7.15) 

\X2 — X3 j;i — X3/ 

Consequently the factor vanishes logarithmically, 

||._L +0(olna) . (7 . 16) 

This again spoils the possible dominance of the zero modes, the result read- 
ing 



1 



a>i 

«2 



Xl - x% 



1 



Xl — £3 

a 3 ,a 4 ~ — | — . (7.17) 
a in a 



So although becoming singular in the limit ( 7.14| ) the zero mode contribu- 
tions will still be sub-dominant in the basis eq. Q6.3| ). The consequences for 
the iterative procedure in the d.s.l. are the same like for the merging cuts, 
where in eqs. ( |7.12j ) and ( |7.13| ) the index 2 has to be replaced by 4. 

An even different kind of continuum limit can be found, when the pro- 
cess of merging and shrinking is put together. In this limit the zero mode 
contributions in the basis will no longer be sub-dominant. Starting with the 
following parametrisation, 

X2 = £3 + av , v > , 
%3 = x% 

X4 = X3 — a^i , n > , 

p = x c 3 + a-K , (7.18) 
the modulus and also the ratio j^ffi will stay finite, 



v 



k 2 = — + 0(a) . (7.19) 
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The modulus may thus reach any value in (0,1). For the special choice v = fx, 
it holds that k 2 = k' 2 = i. The elliptic integrals can then be expressed by 
the gamma function. Using Legendre's relation, which then reads 

2E(k)K(k) - K{k) 2 = | , (7.20) 

and the special value K{-±=) = j^{^{\)) 2 , the ratio of the integrals is given 
by 

The Oj, i = 2, 3, 4 , will therefore scale in the desired way, 

1 



Oil 



xi 



1 / E{k)\ 

«2 = 1 



av \ K{k)J ' 

1 / v + LiE{k) 
a 3 = — -1 + 



av V ix K(k) 

«4 = ^TTT , « = — : — , 7-22 

aixK[k) v + n 



leaving the form of eq. (|6.3| ) unchanged for i = 2, 3, 4. In this kind of d.s.l. 
the iteration merely simplifies, in contrast to the generic case. The only 
terms to be left out in eq. ( 5.2CQ are those with an index 1. 

This ends the short survey of all possible continuum limits to be per- 
formed in the two-cut solution. 



8 Conclusion and outlook 

The results presented here complete the perturbative analysis in of the 
hermitian one-matrix model. All possible solutions classified by the number 
of cuts are derived in an iterative scheme for higher genus contributions 
to the loop correlators, which generalises the results of || for the one-cut 
solution. Explicit expressions were given for the one-loop correlator of genus 
zero and one and for the planar two-loop correlator at different arguments. 
The latter is universal for each number of cuts s, providing a whole set of 
new universality classes. 

The two-cut solution was displayed in full detail, including the genus one 
contribution to the free energy. When investigating the double-scaling limit, 
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the two-cut solution was proven to be equivalent to the one-cut solution in 
the generic case. However, when taking the continuum limit at merging or 
shrinking cuts a different behaviour was revealed. A candidate for this new 
continuum behaviour is still missing. 

Up to now in most investigations of matrix models there has been as- 
sumed a one-cut structure of the corresponding quantities. The results given 
here should make it possible to deal with multi-cut structures also in other 
models, where the loop equation techniques have been applied successfully. 
Just to name a few examples like the 0(n)-model [O] or the supereigenvalue 
model [ff^j, particularly in the complex matrix mod el pCfl most explicit re- 
sults should be accessible, completing the analysis of The hope is beside 
finding new universality classes and continuum limits there as well the ap- 
pearance of new types of critical behaviour, which might still be hidden in 
the full space of solutions. 

Acknowledgements: I would like to thank the Niels Bohr Institute for 
financial support and its warm hospitality, where part of this work was 
being done. In particular I wish to thank J. Ambjorn, C. Kristjansen and 
Yu. Makeenko for many valuable hints and discussions. 
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A jpO) of the boundary conditions 

While taking the derivative ^y(p) of the first type of boundary conditions 
eq. (|4.9| ) is straight forward, the derivative of the second type eq. ( 4.10| ) 

X 2 j 



0= f dXM(X)^Y[^ =i (X- Xi ) , j = l,...,s-l, (A.l) 

is more sophisticated. When applying ^f(p) to M(A) in the integrand, the 
resulting function will no longer be analytic, having an extra pole at p. 
Therefore two things have to be taken into account. First, the contour in 
the expression for M(A) eq. ( |4.6| ) cannot be deformed to infinity but in the 
end of the calculation^], taking it to enclose the poles of the integrand only, 

M(A)=f £ (A.2) 

JCuCx 2m i, , w /rr2s 



Second, when applying J^(p) to eq. ( |A.1[ ) for fixed j, the argument p must 
be excluded to be in the interval of integration, p ^ [%2j+i, %2j]- Otherwise 
the real integral between the cuts does not exist, due to eq. (|5.12|) . The final 
result will have no poles between any of the cuts, allowing for an analytic 
continuation of to p also being between the cuts. 

Taking this into account the derivative -A^ip) in the form of eq. ( |5.11| ) to 
eq. (|Al|) yields 



x 2j 

{ dxi 
J Jc 



duo I 1 



x 2 j+l 



JcuCx 2m \ (p- X) 2 (uj - X) 



2^ dV [P) to-X \uj-xi X-Xi) ) 0(°)(A) 



where now the contours have been deformed back to infinity. The first term 



Otherwise the extra pole at p would be overlooked. 
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can be further evaluated, using partial integration, 

E2 i m-\ , - X2 i 



dp 7 p-A0(°)(A) 2 7 p- X^\X-Xi p-XiJ 0(0) (A) 



which leads to the final form of eq. ( |5.14 ) . 



B The planar two-loop correlator 

The planar two- loop correlator Wo(p,q) with s cuts is derived by apply- 
ing 3p(p) from eq. ( |5.11 ) to the result for Wo(p) eq. (|4.7|) , 



„, , , 11 </> (t %) 11 1 A ,(1), , 

i i , 1 1 y 1 Af. (1)dgi r 

2( 9 -p) 2 4 0(°)(p) £[P-Xi i dV 1 
l^ (0) (g) / 2 1 " 1 | ^g <|t «A 



where the result for the ^(g) eq. ( 5.15| ) has been used. Although this is 



in principle the solution for arbitrary s, it would be nice to have it in a 
form, where the analyticity properties, which follow from its definition, are 
obviously fulfilled. Namely from eq. (|2 .6[) 



W (p,q) = ^(p)W (q) = ±{p)±{q)F , (B.2) 

it is clear, that it must satisfy 

• symmetry W (p,q) = W (q,p) 

• analyticity like Wo (p) ~ (p) , no higher poles 

• asymptotic lim^oo W (p, q) ~ O(^) 
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• regularity lim q ^ p W (p, q) = W Q (p,p) 

i-w = lim e _> (W (p + ie, q)) + W (p - ie, q)) 



(p-g) 2 

The last property can be derived from taking the derivative grf(p) of eq. 
( [4. 3| ), which may be rephrased as 

= 2 Re(W (<?)) , g £ a . (B.3) 

This last point is clearly satisfied by eq. ( [B.ip , whereas the others are not 
obvious. In the case of the two-cut solutions however, where the a^o are 
explicitly known, Wq(j), q) can be cast into a form, where all points are seen 
to be true immediately. Inserting eq. ( J6.4| ) into eq. ( |B.1| ) and using the 
identity 

4 

= 2 Y[(P - x i) + (Q ~P) 2 ({P - x 2 )(p - x 3 ) + (p- xx)(p - x 4 )) 
i=i 

4 

+ (<7-p)$in(P~ x i) ~ (p-»l)(p-^4)(9-aJ2)(?-a:3) 

i=l j^i 

- (q - x\)(q - x 4 )(p - x 2 )(p - x 3 ) (B.4) 
leads to the planar two- loop correlator with two cuts, 

W (p,q) 



4(p 



+ 



l(p- 


x\)(p - 


x 4 )(q - 


x 2 )(q - 


x 3 ) 


(p - 


x 2 )(p - 


X3)(q - 


xi)(q - 


x 4 ) 


Kp- 


- x 2 )(p - 


- x 3 )(q ~ 


- xi)(q - 


~ Xi) 



(p - x\)(p - x A )(q - x 2 )(q - x 3 ) 



1 1 E ^ (I \\ 1 



with 



_ X 2 X 3 X 4 (x 2 - x 4 ) f xix 3 x 4 (xi - x 3 ) 



(xi - X 4 )(x 2 - Xl) (Xl - X 2 )(x 2 - X 3 ) 

+ xix 2 x 4 (x 4 - x 2 ) + xix 2 x 3 (xi - x 3 ) 6 ^ 
(x 2 - x 3 )(x 3 - x 4 ) (xi - x 4 )(x 3 - x 4 ) 

The same can be done for coinciding arguments, using eq. (|5.20| ), which 
reads 

1 4 1 1 4 1 -1 2s 

1 x 1 1 x - 1 1 x Cdj 

W (p,p) = - ^ ^-^5 - g 1- (p _ ^ _ x .) + 4 1. — • ( R7 ) 
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Taking the solution for the at eq. ( ]6.4| ) it can be shown, that 



1 on If 1 1 



^~[P-Xi 4\(p-ii)(p-i4) (p-ar2)(p-aJ3) 

which leads to the final form of eq. ([0]) . It can be checked now, that taking 
the limit q — > p of eq. (|B.5| ) will give back the regular expression Wa(p,p) 
as it should. 



C Elliptic integrals 



The elliptic integrals defined in eq.(|6.2j) can be expressed by elementary 
functions and three fundamental integrals, namely the complete elliptic in- 
tegrals of the first, second and third kind 



K(k) 



dt 



1 



^/(l-t 2 )(l-k 2 t 2 ) 



, E(k) 



dt\ 



'l-k 2 t 2 
1-t 2 



1 

U(a 2 ,k) = J dt— 



(1 - a 2 t 2 )^/(l-t 2 ){l-k 2 t 2 ) ' 



(C.l) 



With the help of the tables in [17| the K{, % = 1, . . . , 4 , are then given by 



X3 



X 



X-2 



Ko = 



/(x 2 -A)(A-x 3 )(A 


- x 4 ) 


/ x\ — A 




a 6 + 2a 4 - 2a 4 *; 2 + 


a 2 *; 2 ) 


- Aa 6 k 2 + 6a 4 k 2 - 


4a 2 *; 2 


/(x!-A)(A-x 3 )(A 


- x 4 ) 


/ x 2 - A 



■L-.i 



X 



-a 6 - 2a 4 + 2a 4 k 2 + a 2 k 2 )E{k) + (-a 4 - 2a 2 + 4a 2 *: 2 - k 2 ) 
(k 2 - a 2 )K(k) + (a 8 - 4a 6 + 6a 4 k 2 - 4a 2 k 4 + & 4 )n(a 2 , fc)l , 
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K A 



%3 

X 

+ ( 

J-2 



(xi - X)(x 2 - A) (A - x 4 ) 
A - x 3 



(3a D - 2a 4 - 2a 4 k 2 + a z k z )E{k) + (k z - a z )(3a 4 - 2a z - k z )K{k) 



.2 ,2 



3a s + 4a b + 4a b k z - 6a A k z + k 4 )U(a z , k) 



(x 1 - X)(x 2 - A)(A - x 3 ) 



X'3 



.(fc 2 



A — X4 

-a 6 + 2a 4 + 2a 4 k 2 - 3a 2 k 2 )E(k) + (-a 4 + 2a 2 
- a 2 )K(k) + (a 8 - 6a 4 k 2 + 4a 2 /c 2 + 4a 2 k A - 3£: 4 )n(a 2 , k) 



Aa 2 k 2 + 3k 2 ) 



(C.2) 



where 



X = 



k 2 = 



3 7 
1 (X 1 - X 3 )2(X 2 - X 4 )2 

4 [x 2 - x 3 ) 2 {x 3 - x 4 ) 
(xi - xa)(x 2 - x 3 ) 
(xi - x 3 )(x 2 - x 4 ) 



a 



x 2 - x 3 
x 2 - x 4 



(C.3) 



The composition into K(k), E(k) and Il(a 2 ,k) is not unique because of 
possible transformations of the modulus k and the parameter a. Using the 
explicit form of eq. (C.2), the coefficients a.% in the solution for thejfe(p) 
eq. (|6.3j ) can be shown to have the simple form, that is given in eq. fl6.4|) . 
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